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Abstract 

We study the production of A hyperons in p + A reactions on the basis of a 
BUU transport approach from 1.1 to 1.9 GeV and evaluate the properties of the 
hypernuclei produced in particular with respect to their momentum distribution 
in the laboratory frame. Due to elastic AA^ scattering large cross sections for 
the production of heavy hypernuclei in the order of 100 - 400 jib are predicted 
for p + Uatl.5-1.9 GeV laboratory energy. Whereas the K~^Y production 
channels are expected to be only slightly modified in the nuclear medium, the 
antikaon production should be enhanced substantially due to large attractive K~ 
selfenergies in dense matter. We predict an enhancement of the inclusive K~ yield 
in p + 208p]^ collisions of a factor of ~ 10 at 2 GeV laboratory energy. 



1 Introduction 

The investigation of strangeness in hadrons and nuclei has become an exciting and chal- 
lenging field of research in the last years IQ. Here, the open questions reach from the 
spectroscopy of hypernuclei to the formation of strange hadronic matter in neutron stars 
or ultrarelativistic nucleus-nucleus collisions. Whereas earlier studies were devoted to 
the production aspects of hyperons and strange mesons, nowadays, the interest is moving 
towards the properties of strange particles in a nuclear environment. 

A-hypernuclei are especially well suited for in-medium strangeness investigations 
since the A-hyperon has a long lifetime as compared to the nuclear scale. The 'free' 
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A-hyperon decay is purely mesonic (A — > Nn) whereas in the nuclear medium also non- 
mesonic decay channels {AN — >■ NN) are possible. The competition between these decay 
modes is expected to provide information about the hypernuclear structure (e.g. the A 
and pion selfenergies in the nuclear medium) 0. Since the lifetime of the hypernucleus 
depends on the corresponding widths of these decays , i.e. t = h / [ Tmes + ^nonmes], 
it is of considerable interest to measure the lifetime of hypernuclei as a function of their 
mass, particularly for heavy systems, where Tmes is strongly suppressed due to Pauli 
blocking such that Tnonmes is directly related to r~^. 

This particular question has been investigated with antiproton beams on ^°^Bi and 
^^^U targets in However, the lifetimes of ^°^Bi and ^^^U from this experiment are 
quite different contrary to simple phase-space expectations. Alternatively, one might 
investigate the production of heavy hypernuclei by means of proton + nucleus reactions. 
In fact, the recent studies on the (p , K+ ) reaction confirm a quite substantial production 
of associated A-hyperons Q leading to production cross sections for A-hypernuclei in 
the order of a few 100 fxb for p + Pb at 1.5 - 1.9 GeV |, j^. 

Apart from the lifetime of the A hyperon in the medium also the kaon and antikaon 
properties should change due to interactions with the nuclear environment |P, |^, 0. 
Whereas the potential is expected to be slightly repulsive in the medium, the K~ 
meson should see a sizeable attractive potential in dense matter such that their produc- 
tion should be enhanced at finite baryon density. In case of heavy-ion collisions a first 
exploratory study has been performed by Li, Ko and Fang ||TT| with the result that large 



attractive K~ potentials are needed to explain the experimental spectra from ||T2[ for 
Ni -|- Ni at 1.85 AGeV. A more systematic transport analysis - including all possible 
production channels - on this question has been performed in Refs. [|l^, 0] in compari- 



son to the recent data from the KaoS Collaboration |T^. In fact, the experimental data 
indicate an attractive K~ potential of about -100 MeV at normal nuclear matter density 
(po ~ 0.16 fm~^) whereas the potential at density po is in the range between and 
+30 MeV |jl3|, 0]. Such effects should also be seen and probed experimentally in p + A 



reactions at subthreshold energies. 

In continuation of our earlier work 0, |^ we here present a transport (BUU) analysis 
of proton + nucleus collision events, including the production channels pN —>■ NYK^ 
and ttN YK^ (Y = A, S), while nonperturbatively taking into account the rescat- 
tering of the hyperons with nucleons and the hyperon propagation in the nuclear mean 
field. 



2 Ingredients of the transport approach 

In this contribution we perform our analysis along the line of the HSD^ approach |]l7 



which is based on a coupled set of covariant transport equations for the phase-space 
distributions fh{x,p) of hadron h, i.e. 

{(n, - u^df^u;: - Midiul) a^f + (n.a^t/,^ + Midiu^) d>;] Mx,p) 



^Hadron String Dynamics 
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= J d2d3dA...[G^G]u^u..ATl + U2-U3-U^...) 

h2h3h4,... 

X {fh,{x,P3)fhAx,P4)fh{x,p)fh2{x,P2) 

- fh{x,p)fh2{x,P2)fh3{x,P3)fh^{x,p4)}... . (1) 

In Eq. (|l]) U^{x,p) and Uf^{x,p) denote the real part of the scalar and vector hadron 
selfenergies, respectively, while [G~^ G]i2-,34.. 6^(11 + 112 — lis — 114 . . .) is the 'transition 
rate' for the process 1 + 2-^3 + 4 + ... which is taken to be on-shell in the semiclassical 
limit adopted. The hadron quasi-particle properties in (P are defined via the mass-shell 
constraint 

^(n^n'^-Mf) , (2) 

with effective masses and momenta (in local Thomas-Fermi approximation) given by 

M*ix,p) = M,, + U^{x,p) 

U>^{x,p) = p^-U\:{x,p) , (3) 

while the phase-space factors 

fh{x,p) = l±fh{x,p) (4) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depend- 
ing on the type of hadron in the final/initial channel. The dots in Eq. (|1|) stand 
for further contributions to the collision term with more than two hadrons in the fi- 
nal/initial channels. The transport approach ^ is fully specified by U^{x,p) and 
Uh{x,p) {fi = 0,1,2,3), which determine the mean-field propagation of the hadrons, 
and by the transition rates . .) in the collision term, that describes the scatter- 

ing and hadron production/absorption rates. 

The scalar and vector mean fields and ?7/^ for nucleons are taken from Ref. 



the hyperon mean fields are assumed to 2/3 of the nucleon potentials. In the present 
approach, apart from nucleons, A,iV(1440), A^(1535), A, S with their isospin degrees 
of freedom, we propagate explicitly pions, K~^,K~, and r^'s and assume that the pions 
as Goldstone bosons do not change their properties in the medium; we also discard 
selfenergies for the ?7-mesons which have a minor effect on the K~^,K~ dynamics. The 
kaon and antikaon potentials, however, have to be specified more explicitly. 



2.1 K^jK selfenergies 

There are a couple of models for the kaon and antikaon selfenergies ||^, 0, which 
differ in the actual magnitude of the selfenergies, however, agree on the relative signs 
for kaons and antikaons. Thus in line with the kaon-nucleon scattering amplitude the 
potential should be slightly repulsive at finite baryon density whereas the antikaon 
should see an attractive potential in the nuclear medium. Without going into a detailed 
discussion of the various models we adopt the more practical point of view, that the 
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actual and K selfenergies are unknown and as a guide for our analysis use a linear 
extrapolation of the form, 

m*^{pB)=mll^ + a^^, (5) 

with a ~ -0.2 for antikaons and a ~ 0.06 for kaons (or a = for the bare kaon). Our 
choice {a ~ -0.2) leads to a fairly reasonable reproduction of the antikaon mass from 
Ref. p (thin solid line in Fig. 1) and the recent results from Waas, Kaiser and Weise 
p!0| (thick solid line in Fig. 1). We note that the dropping of the antikaon mass is asso- 
ciated with a corresponding scalar energy density in the baryon/meson Lagrangian, such 
that the total energy-momentum is conserved during the heavy- ion collision (cf. ]T7|). 
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Fig. 1: The , K~ mass as a function of the baryon density in units of po ~ 0.16/m~^ 
according to Kaplan and Nelson |^ (thin solid lines) and Waas, Kaiser and Weise [10 
(thick solid lines). 



2.2 Perturbative treatment of strangeness production 

The calculation of 'subthreshold' particle production is described in detail in Ref. |]TB| 
and has to be treated perturbatively in the energy regime of interest here due to the 
small cross sections involved. Since we work within the parallel ensemble algorithm, each 
parallel run of the transport calculation can be considered approximately as an individual 
reaction event, where binary reactions in the entrance channel at given invariant energy 
^/s lead to final states with 2 (e.g. K^Y in tiB channels), 3 (e.g. for K^YN channels 
in BB collisions) or 4 particles (e.g. KKNN in BB collisions) with a relative weight 
Wi for each event i which is defined by the ratio of the production cross section to 
the total hadron-hadron cross section^. We thus dynamically gate on all events where 
a K^Y or K^K^ pair is produced initially. Each strange hadron then is represented 
by a testparticle with weight Wi and propagated according to the Hamilton equations 

^The actual final states are chosen by Monte Carlo according to the 2, 3, or 4-body phase space. 
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of motion. Elastic and inelastic reactions with pions, ?7's or nonstrange baryons are 
computed in the standard way and the final cross section is obtained by multiplying 
each testparticle with its weight Wi. In this way one achieves a realistic simulation of the 
strangeness production, propagation and reabsorption during the heavy-ion collision. 

For further details of the model and the explicit parametrizations of the various 
reaction cross sections we refer the reader to Refs. 0, O, O, 



14 1, respectively. 



3 Recoil momentum distributions 

Since in the transport calculations the four-momenta of all hadrons are propagated in 
time we are capable to extract for events, leading to hypernucleus formation, the average 
properties of the residual hypernucleus. For this purpose we compute as a function of 
time those particles (essentially nucleons) that have left the residual heavy fragment at 
position R, i.e. 

|r, - R| > i?^ + 2/m, (6) 

1 /3 

where Ra = 1.2 fm denotes the radius of the target with mass number At. Now 
let the number of particles emitted be Np{t). For each parallel ensemble we can then 
evaluate the fragment's average mass number, its excitation energy, three-momentum 
and angular momentum by exploring the conservation of total energy, mass number, 
momentum and angular momentum: 

<E*> (t) = Et,t - J2 ^Jp] + Mf - Mres - E,oui -M^ + Mn- Ek, 
<AF>it) = At + l-Np{t), 

< p > (t) = Ptot - Pj 

7Vp(i) 

< L > (t) = Ltot - E ^^jit) X Mt)- (7) 

i=i 

In eq. (|^ Mres denotes the mass of the 'residual nucleus', Ecoui stands for the Coulomb 
energy between the emitted particles and the 'residual nucleus', while Ek represents the 
total energy of the i^'"'"-meson. 

All quantities in depend explicitly on time t due to the continuous evaporation 
of particles from the final compound system. Since we will follow the further decay 
chains by statistical model codes, the actual transition time for the connection of the 
BUU and the statistical model calculation is of no significance as long as the system has 
left the nonequilibrium phase of the reaction and achieved statistical equilibrium. We 
have checked that it is sufficient to trace the history of each ensemble of events within 
BUU up to 150 fm/c §. 

As an example for the momentum distribution of hypernuclei we show in Fig. 2 
the computed distributions in beam direction (dN/pz) and perpendicular to the beam 
(dN/dpj.) for p + ^^^f/ at Tiab = 1-5 GeV from the BUU calculation (dashed lines) as 
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well as after particle evaporation (solid lines) computed via PACE2 [jT8|, |^ . We note that 
the particle evaporation only broadens the momentum distribution somewhat, but does 
not change the average momentum per particle in the beam direction. 
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Fig. 2: Longitudinal momentum (pz/A) and transverse momentum (px/A) distribution 
of A-hypernuclei at Tiab = 1.5 GeV for p + "^^^U with evaporation (solid histograms) and 
without evaporation (dashed histograms). 



Since the computed momentum distributions enter as an important ingredient in 
the experimental analysis ||19| , |20| based on the recoil shadow method, we have to investi- 
gate the accuracy of the transport calculations with respect to the momentum transfer in 
proton-nucleus reactions. In this respect we show in Fig. 3 the longitudinal momentum 
distribution of the residual nuclei from the BUU calculation (solid histograms) - without 
gating on hypernuclei - for p + ^^^V at Tiab = 475 MeV, 1.0 GeV, 1.5 GeV, and 2.9 GeV 



in comparison to the data of Fraenkel et al. |^ (full squares) and Kotov et al. [^. The 
good agreement with the data in this wide kinematical regime demonstrates the relative 
accuracy of the transport approach which should be of the same quality when gating on 
events with hypernucleus formation (cf. Fig. 2). 
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Fig. 3: Comparison of tlie longitudinal momentum distribution for p + "^^^U at Tiat = 
475 MeV (Ihs. top), 1.0 GeV, 1.5 GeV, and 2.9 GeV from the BUU calculation (solid 
histograms) with the data from Fraenkel et al. |2l| and Kotov et al. (full squares). 
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4 and K production in p + A collisions 



The inclusive production of kaons and hypernuclei in p + A collisions has been studied at 
subthreshold energies to a large extend in Ref. |^ and does not have to be reviewed 
here. The novel aspects are those related to antikaons and their selfenergy or potential 
inside heavy nuclei. Since antikaon production has a threshold of about 2.5 GeV in free 
nucleon-nucleon collisions, we have to address the question of kaon production at Tiab ~ 
2 - 3 GeV, too. 

In this respect we first compare our calculations (solid histograms) with the ex- 
perimental spectra for p + Pb at 2.1 GeV from Schnetzer et al. in Fig. 4 
without including any medium modification of the kaons (in line with the hypothesis 
of Kaplan and Nelson (cf. Fig. 1)). Indeed, the experimental spectra are described 
quite accurately as in case of nucleus-nucleus collisions at SIS energies without any 



selfenergies such that a measurement oi K / ratios at the same bombarding energy 
becomes interesting due to its sensitivity to the antikaon potentials. 
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Fig. 4: The calculated spectra for p -|- Pb collisions (solid lines) at 2.1 GeV in 
comparison to the data from Schnetzer et al. at various angles Oiah in the laboratory. 
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Employing the various production channels for antikaons {pN NNKK, irN 
NKK, ttY NK, YN NNK) as well as antikaon absorption {KN ttY) from 
Ref. [|13| we have performed detailed calculations on the and K~ spectra in p + 
^°^Pb collisions from 2-3 GeV bombarding energy. The first predictions for the total 
K^/K^ ratio are presented in Fig. 5 (integrated over all momenta) as a function of 
Tiab for a bare antikaon mass (open squares) and an in-medium antikaon mass (full dots) 
according to Eq. (5) with a = -0.2 which corresponds to an attractive antikaon potential 
of about -100 MeV at po- We find the K~ / ratio to increase from 2x10""^ to 7xl0~^ 
in this energy regime for the bare antikaon case; however, when including the attractive 
K~ potential the ratio is enhanced by about a factor of 10 at Tiab = 2 GeV and by a 
factor of ~ 2 at Tiab = 3 GeV. 




2.0 2.5 3.0 
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Fig. 5: The calculated K~ / ratio for p + ^osp]^ reactions from 2 to 3 GeV; with an 
attractive antikaon potential (full dots); without antikaon potential (open squares). 



5 Summary and conclusion 

In this study we have presented a BUU transport analysis of K^Y and K^K^ production 
in proton + nucleus collisions at COSY energies employing the elementary production 
processes from our earlier work [T^ . Due to a large fraction of A hyperons from the 
secondary process vriV — »• K^K, which leads to hyperons with moderate momenta in the 
laboratory system, the p-|-nucleus reaction efficiently produces heavy hypernuclei and 
cross sections of about a few 100 iib are expected for p + ^osp]^ ^.^ g _ j^ g Qgy g]^ 
which is in line with present experimental data so far |l^, |2^ . 

We have, furthermore, shown that the transport approach also reliably describes 
the momentum transfer to the target nucleus in a wide kinematical regime from 475 
MeV to 2.9 GeV for heavy nuclei such that the computed momentum distributions for 
hypernuclei, that are needed for the experimental analysis |1T9| , pO| based on the recoil 
shadow method, are expected to have the same accuracy. 
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In addition, we have explored the possibihty to measure the K~ potential in finite 
nuclei via their production cross section relative to mesons. In fact, p + ^osp]^ 
reactions at 2.0 GeV indicate an increase of the K'^ cross section about a factor of 10 
when including a potential of about -100 MeV at po (cf. Fig. 1) in line with the 
chiral Lagrangians of |TD[ and the analysis in Ref. |T^. This enhancement should be 
clearly seen in the next generation of experiments. 
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